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1 Introduction 

The multiparticle spectra of secondaries produced in high energy collision processes 
are the most abundant sources of our knowledge of the dynamics of such processes. 
Among others, two features emerging from the analysis of these spectra are of partic- 
ular interest: (i) the so called intermittent behaviour observed in many experiments 
in the analysis of factorial moments of spectra of produced secondaries and (ii) the 
Bose-Einstein correlations (BEC) observed between identical particles. Whereas the 
former seems to indicate the existence of some (multi) fractal structure of the pro- 
duction process |1] the latter are established as, by now the most important source 
of our knowledge on the space-time aspects of the multiparticle production processes 

§■ 

Some time ago it was argued || |] that, in order to make both effects compati- 
ble with each other, the emitting source should fluctuate in size in a scale-invariant 
(i.e., power-like) way. This can be achieved in two ways: (i) either the shape of the 
interaction region is regular but its size fluctuates from event to event according to 
some power-like scaling law (ii) or the interaction region itself is a self-similar fractal 
extending over a very large volume || Q. 

In this work we would like to investigate in more detail to what extent the BEC is 
sensitive to the possible space-time fractality of the emission source. To this end we 
shall use a simple self-similar cascade process in which the final particles are pro- 
duced in the sequential two-body decays of some original mass M. For our purpose 
we shall extend it by introducing the simple (classical) space-time development of 
the cascade and by adding the kind of BEC "afterburner" along the lines advocated 
recently in [J/j. 

It is widely expected that every cascade model has automatically built in the 
intermittent behaviour of spectra of observed particles ||. Although this statement 
is true and obvious for the models based on random multiplicative processes in some 
chosen observed variables (like energy, rapidity or azimuthal angle) it is highly non 
trivial in the case of cascades (or multiplicative processes) proceeding in variable (s) 
not directly mesurable but nevertheless of great dynamical importance (as, for ex- 
ample, masses of some intermediate objects occuring during the production process 
HI). In the purely mathematical case, where cascade process proceeds ad infinitum, 
one eventually arrives at some space-time fractal picture of the production process. 
However, both the finite masses of produced secondaries and limited energy (or 



mass M) stored in the emitting source prevent the full and distinct development of 
such a fractal structure ||. One must therefore be satisfied with only some limited 
and mostly indirect presence or signals of such structure. If established it would, 
however, be important for our knowledge of the dynamics of the multiparticle pro- 
duction process. 

Such a fractal structure in phase-space can generate a similar structure in the 
space-time picture of the hadronization process. Our aim here is to demonstrate to 
what extent they influence BEC In the next two Sections we shall then provide, re- 
spectively, the phase-space and space-time characteristics of a simple cascade model 
used for that purpose. Section 4 contains our main results showing the BEC fea- 
tures emerging from our model. Section 4 contains a summary of our results and 
conclusions. 



2 Phase-space characteristic of the cascade model 
used 

We shall model the emitting source of mass M by the usual (1 — > 2) random cascade 
process employed already in ||, M — ► Mi + M 2 , in which the initial mass M 
"decays" into two masses M lj2 = fc li2 ■ M in such a way that k\ + k 2 < 1, i.e., a 
part of M equal to (1 — k± — k 2 )M is transformed to kinetic energies of the decay 
products Mi t 2- The process repeats itself (see Fig. 1) until M\ )2 > // (/i being the 
mass of the produced particles) with successive branchings occurring sequentially 
and independently from each other, and with different values of fc lj2 at each branch- 
ing, but with energy-momentum conservation imposed at each step. For different 
choices of dimensionality D of our cascade process (provided by the restrictions for 
the possible directions of flights of the decay products in each vertex) be it D = 1 or 
D = 3 dimensional (isotropic) and for different choices of the decay parameters k\$ 
at each vertex, we are essentially covering an enormously vast variety of different 
possible production schemes ranging from the essentially one-dimensional strings to 
thermal-like fireballs. 

Of special interest is the case of a one-dimensional cascade for which one can 
provide analytic formulae for the rapidities Y 12 of the decay product at each vertex 
given in the rest frame of the parent mass in this vertex. They depend solely on the 
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decay parameters at this vertex, fc 12 : 
Fi = ± In 
Y 2 = Tin 
where 

Two limiting cases can be distinguished here: (i) - totally symmetric and (ii) 
maximaly asymmetric cascades. In the case of a totally symmetric cascade decay 
parameters are equal and the same for all vertices, fc 12 = k. In this case the finally 
produced particles occur only at the very end of the cascade process and the amount 
of energy allocated to the production is maximal. Because the number of possible 
branchings characterizing the length of the cascade is equal to L max = In — / In £ 



(where \x = y m 2 , + (pr) 2 ), the multiplicity of produced secondaries is given by the 
following formula: 

r (M\ dF In 2 



According to || 10 1 the exponent dp is formaly nothing but a generalized (fractal) 



dimension of the fractal structure of phase space formed by our cascade. The utility 
of such notion is, however, greatly reduced because of the necessary limited length 
of our cascades . Notice a kind of scaling in ([|) where, for a fixed ratio — the ob- 
served multiplicity N s depends solely on the decay parameter k. The characteristic 
power-like behaviour of N S (M) in (0) is normally atributed to thermal models. For 
example, for k — \ one has N s ~ Ma, which in thermal models would correspond 
to the ideal gas equation of state with velocity of sound °o = ~~m EUl- The same 
behaviour (on average) is obtained for fc 12 chosen randomly from a triangle distri- 
bution P(k) = (1 — k) a with a ~ 1 which will be therefore used in all our numerical 
calculations. The fact that iV s E (2, —J means that decay parameters are limited 

tofce(£,i). 

For maximally asymmetric cascades k\ = ■£ and k 2 = k, i.e., at each step one 
has always a single final particle of transverse mass \x produced against some recoil 
mass Mi = kM: M — » \x + M\. The amount of kinetic energy allocated to the 
produced secondaries is now maximal. The corresponding rapidities Y\ and Y 2 are 
now given by eqs. ([I]) with, respectively, k\ — jr and k 2 = k. Also here the resultant 
multiplicity iV a is given by the length L max of the cascade, which is limited by the 
condition M/ = k l M > //. It means therefore that L max = In— /ln^ and the 



corresponding multiplicity is 

1 M 

N a = 1 + L max = 1 + —3- • In — . (3) 

lni n 

The energy dependence of N a is now logarithmic (in thermal models it would cor- 
respond to a one dimensional fireball with cq — > 1 ]12[) but the same kind of scaling 
as in eq. (|2[) is present also here. 

Our simple model therefore covers all possible energy dependencies of the mul- 
tiplicities of produced particles which depend solely on decay parameters ki^ of the 
cascade and scale in the ratio — . This remains true for both one and three dimen- 
sional cascades. 

In Fig. 2 one can see examples of rapidity distributions j^^f- calculated for sym- 
metric and asymmetric D — 1 cascades discussed above. They are compared there 
with most probable distributions in one-dimension obtained by means of information 
theory arguments [1^ 



frr(v) = -z ■ exp [-(3 ■ fi coshy] , (4) 

where Jdyf(y) = 1 (what defines Z) and j3 = /3(M,N) is the corresponding la- 
grange multiplier ensuring proper conservation of energy-momentum in the case 
when N particles, each of transverse mass /i, are produced from the source of mass 
M. Contrary to the case of production via the cascade process, nothing is now 
said about the production mechanism. It is just tacidly assumed that all produced 
particles occur, in a sense, instantenously in the whole allowed phase space with 
the weights provided by eq. (|4]), which was obtained by the maximalization of the 
suitably defined information entropy corresponding to the production process under 
consideration [[y|. It occurs that for a wide range of M and iV the quantity J3 = /3|l 



is (almost) constant as a function of energy per particle m = j^, i.e., also here one 
encounters a kind of scaling, namely that -^ — ~ F ( z 
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The shape of the multiplicity distribution, P(N), in our case of a source with 
fixed ratio — is given by distribution P(kip) of decay parameters kip- We shall use 
a simple triangle form for it (as already mentioned above) P(k) = (l — k) a , which for 
a ~ 1 provides the commonly accepted energy behaviour of the mean multiplicities 
N(M) ~ j\^f0.4-o.5 ag di scussec i above. The example of P{N) for D = 1 cascades are 
shown in Fig. 3 (P(N) for D = 3 cascades are the same). They exemplify three dif- 
ferent choices of the ratio — (^ = 33.3, || = 133.3 and ^ = 333.3, respectively). 



In Fig. 4, we show the behaviour of scaled ("horizontal" [ll]) factorial moments 
Fi (for I = 2,3) calculated for a one-dimensional cascade in rapidity space as a 
function of number of bins nun = Y/Sy (where Y is taken as the corresponding 
rapidity range for the corresponding mass M and Sy denotes the bin size considered) : 

-i)---(m-l + i)\ 




(Ny 



(5) 



The N produced particles are distributed among nu n bins with n* particles in the i th 
bin, i.e., Y%=i n i = N. The average is over all events. It is interesting to note that 
our results, although obtained for essentially the same type of cascade as discussed 
in 0, apparently demonstrate much stronger intermittency signal than the experi- 
mental one shown there. However, no fit to the data was attempted in our case as 
we are concerned with the properties of a single elementary source only (leaving the 
problem of their distribution in mass P(— ) aside). On the other hand, moments F\ 
in [|] were in reality not calculated but deduced from experimental data by means of 
some simple formula obtained from general (mathematical) fractal analysis of sym- 
metric cascade processes. The aim was to deduce from them the fractal dimensions 
dp = In 2/ In r of cascade process considered. As a result the corresponding decay 
parameter in || turns out to be very large, k ~ 0.45, leading to (N(M)) ~ ^0.8^0.9 
instead of expected (N(M)) ~ M 0A ^ - 5 as discussed above IT 



3 Space-time characteristic of the cascade model 
used 

We shall now endow our cascade in phase-space with space-time elements (not ad- 
dressed in ||). To this aim we introduce some fictitious finite "life time" t for each 
vertex mass Mi, which is allowed to fluctuate according to some prescribed distribu- 
tion law T(t). This procedure is a purely classical one, i.e., we are not treating Mi as 
resonances, as was done, for example, in [plj on another occasion. Instead, they are 



regarded to be real particles with masses given by the corresponding values of decay 
parameters k\$ and with the respective velocities equal to (3 — -^ ((£1,2; -Pi, 2) are 
the energy-momenta of the corresponding decay product given in the rest frame of 
the parent mass in each vertex). The energy- momentum and charges are strictly 
conserved in each vertex separately (this is another difference with the information 
theory approach [13| where such conservation laws are imposed on the whole process 
instead). As for a decay/branching law we shall choose it in the simplest possible 
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exponential form: 

T(t) = --exp — (6) 

T L tJ 

It is straightforward to get our cascade model in the form of a Monte Carlo 
code. The main features of a one-dimensional case has already been demonstrated 
above. The only difference between one and three-dimensional cascades is in the 
fact that, whereas in the former decay products can flow only along one, chosen 
direction, in the latter in each vertex the flow direction is chosen randomly from 
the isotropic angular distribution. To allow for some nonzero transverse momentum 
in the one-dimensional case we are using the transverse mass \x = 0.3 GeV. For 
the three-dimensional cascade /i is instead set simply to the pion mass, /i = 0.14 
GeV. In every case all decays are described in the rest frame of the corresponding 
parent mass in a given vertex. To get the final distributions, one has to perform a 
number of Lorentz transformations to the rest frame of the initial source mass M. 
As output we are getting in each run (event) a number Nj of secondaries of mass 



/i with both defined energy-momenta 


E j = vV + Pi Pi 


and space-time 

» — 1 AT. 


coordinates [i,-; fj] i _ 1 N , of the last branching (i.e., the coordinates of birth of each 


particle). 







Fig. 5 shows densities p(r) of points of production for all cases investigated here: 
for D = 1 and 3 dimensional cascades with both constant and mass dependent 
evolution parameter r and for three choices of the source mass, M = 10, 40 and 
100 GeV. As one can see the widely expected (cf. || ||) power-like behaviour of 
cascading source 

p{r) ~ ( -) r>r Q , (7) 



is seen only (if at all) for r > r , i.e., for radii larger then some (not sharply defined) 
radius r , value of which depends on all parameters present here: mass M of the 
source, dimensionality D and evolution parameter r of the cascade. Below tq the 
p(r) is considerably bended, remaining even almost flat for D = 1 cascades. Only 
for rapidly developing cascades (i.e., for r ~ -g) in D = 3 dimensions, the expected 
scaling sets in almost from the very beginning and it practically does not depend 
on the mass of the source. For the limiting case of M = 100 GeV the correspond- 
ing values of parameter L vary from L = 1.89 and L = 1.86 for r = 0.2 and 1/M 
for one dimensional cascades to L = 2.78 and L = 2.8 for three dimensional cascades. 

The shapes of p(r) scale in the ratio f — ) in the same way as the multiplicity 
distributions P(N) discussed before. As desidered power-like behaviour 0-0 sets 



in (at least approximately) only for long cascades (large values of ( — )) and/or for 
fast ones (small values of r), it remains therefore to be checked whether (and to 
what extent) such conditions are indeed met in the usual hadronic processes. This 
point is, however, outside of the scope of the present work. 



4 Cascades and BEC 

Let us proceed now to our main point, namely to the question of whether one can see 
in BEC some special features which could be attributed solely to the the branchings 
and to their space-time and momentum space structure. At first glance, the answer 
seems to be plainly negative as it is easy to check that the function 



ft(« = |*-i>,|)- dN( *" p ' : 



dN(p t ) dN{pj) 



does not show any structure of BEC type. It is also true if we endow our cascade pro- 
cess with the production of charges of the type: {0} — ► {+} + {—}, {+} — ► {+} + {o} 
and {— } — > {0} + {— }. In this case Ci calculated for like charge pairs also does not 
show any correlations. That is, however, to be expected, because the only way to 
have (|D showing "primordial" BEC is to introduce them from the very beginning, 
for example in the way done recently in ]nj. Using the same information theory 



approach as in [J^j but adding a new piece of information, namely that the produced 
particles are mostly bosons and as such they should be grouped together as much 
as possible in the phase space cells, the authors of [|15|] indeed obtained a substan- 
tial BEC signal, namely Ci > 1 for pi — > pj in (§). No space-time structure was 
discussed in I|l5fl, however. 



We cannot follow this prescription here without invoking a kind of extremely 
difficult to formulate (or calculate) special final state interactions between produced 
secondaries. Our cascade is supposed to mimic the production process in its de- 
velopment, whereas the information theory procedure of |0| makes no statements 
whatsoever about the development of the proces as such. It only provides the least 
biased and at the same time most probable distributions, limited only by imposed 
constraints of the energy-momentum and (mean) number of particles conservation 
(reflection of which are the two lagrange multipliers, (3(M, N) and fi(M, N), repre- 
senting in terminology of the usual thermal models the "inverse temperature" and 
"chemical potential", respectively. We could, in principle, use the pseudopotential 
method as, for example, advocated long ago in [|I7| , but this causes changes in the 
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particle distributions and/or destroys the energy-momentum balance which has to 
be later restored in a more or less ad hoc way and it does not use the information 
on the space-time structure of our results. 

An open question is the possible existence of phase factors in every branch point, 
which would endow each particular branch and through it also the finally produced 
secondaries with some specific, possibly path dependent phases. This question is, 
however, left open here and it is understood that they are all set equal to unity. They 
would be important to BEC correlations influencing especially values of C 2 (Q = 0), 
i.e., the so called degree of coherence/chaoticity A. We shall return to this problem 
elsewhere. 

Because our aim is not data fitting but checking if, and to what extent, the 
BEC via its C 2 observable is sensitive to different choices of the cascade processes 
provided by different sets of parameters, we have decided to use the ideas of the 
BEC "afterburners" advocated recently in J7j. And because we are not so much 
interested in particular values of the "radius" and "coherence" parameters R and 
A, but in the systematics emerging from our study, we shall use for this exploratory 
research the most primitive, classical version of such "afterburner" . The procedure 
we use is therefore very simple. After generating a set of i = 1, . . . , Ni particles for 
the I th event we choose all pairs of the same sign and endow them with the weight 
factors of the form 

C = 1 + cos [(n - rj)(pi - Pj)] (9) 

where r; = (tj,rj) and Pi = (Ei,pi) for a given particle. 

The results obtained from N event = 50000 events are presented in Fig. 6. They 
are displayed for the same sequence of parameters M (mass of the source), D (its 
dimensionality) and r (the evolution parameter) as in Fig. 5. The characteristic 
feature to be noted is a substantial difference between D — 1 and D = 3 dimensional 
cascades both in the widths of the CziQ) and their shapes. Whereas the former are 
more exponential-like the latter are more gaussian-like with a noticeably tendence 
to flattening out at very small values of Q. Also values of intercepts, C 2 (Q = 0), are 
noticeable lower for D = 3 cascades. There is also a difference between "slow" (con- 
stant r) and "fast" (r ~ -g) cascades, especially for D = 1 ones. The former lead 
to substantially different shapes in this case. In D = 3 this effect is not so visible, 
although it is also present. The length of the cascade (i.e., the radius of the produc- 
tion region, cf. discussion of density p before) dictates the width of ^(Q). However, 



the ( — ) scaling observed before in multiplicity distributions and in shapes of source 
functions is lost here. This is because C2 depends on the differences of the momenta 
p — a coshy, which do not scale in — . The flattening mentioned above together 
with 6*2(0) < 2 for D = 3 cascades are the most distinctive signature of the fractal 
structure combined with D = 3 dimensionality of the cascade. The correlations of 
the position-momentum type existing here as in all flow phenomena are, in the case 
of D = 3 cascades, not necessarily vanishing for very small differences in positions 
or momenta between particles under consideration. The reason is that our space- 
time structure of the process can have in D = 3 a kind of "holes", i.e., regions in 
which the number of produced particles is very small. This is perhaps the most char- 
acteristic observation for fractal (i.e., cascade) processes of the type considered here. 



5 Summary and conclusions 

In this work we have addressed the problem of the possible space-time fractal struc- 
ture of the hadronic production process. It is complementary to the possible (multi) 
fractality leading to fluctuations in the multiparticle distributions and to the possible 
fractality claimed to exist already on the level of hadronic structure [pl|. Although 
there is a vast literature concerning the possible (multi)fractality in momentum 
space |19j its space-time aspects are not yet fully recognized with 0-0 remaining 



so far the only representative investigations in this field. Our aim was to extent 
this investigation a bit further by essentially repeating ideas proposed in ||-0 in a 
numerical form that allowed us to check in more detail the conjectures made there, 
showing the limits of their applicability. 

Our simple model posseses all features called for in - : it shows both inter- 
mittency in the phase space (demonstrated in the limiting case of one-dimensional 
cascade explicitely in Fig. 4) and (approximate) power law distribution of the pro- 
duction points in the space-time (cf. Fig. 5). As we have more constraints on the 
phase space behaviour imposed by the, for example, expected energy dependence 
of the multiplicity (N), we have not much freedom in choosing decay parameters 
k. Our distribution P(k) is surely not the only possible one. However, whatever 
shape we choose for P(k) it should reproduce the expected energy dependence of 
multiplicity, (N(M)) ~ ^0.44-0.5 Therefore the only really free parameters in our 
study were the time evolution parameter r and dimensionality of the cascade, D. 
Here only two extreme values of D = 1 and D = 3 were studied and two also, in a 
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sense, extreme behaviours of r = const and r ~ -^ were used. All three: M, .D and 
r were found to influence the C^Q) observable characterising BEC, cf. Fig. 6. 

We therefore conclude that BEC are, indeed, substantially influenced by the fact 
that our process is of the cascade type as was anticipated in @-f|, although proba- 
bly not to the extent expected there (which, however, has not been quantified there). 
However, in practical applications, i.e., in the eventual fitting of experimental data, 
there are many points which need further clarification. The most important is the 
fact that data are usually collected for a range of masses M and among directly 
produced particles are also resonances. Therefore, one has first to specify the form 
of the distribution P{— ) which will influence to some extent our results. In partic- 
ular changes in \i due to the production from resonances will shorten our cascade 
considerably. The possible effect can be to some extent deduced from our results 
by comparing data with M = 100 GeV with those with M = 40 GeV and M = 10 
GeV. The BEC will be effective only in conjunction with precise studies of distri- 
butions in the phase-space (like P(— ), intermittency and momentum and rapidity 
distributions. These studies should to some extent fix the distribution of decay pa- 
rameters P(kij). Only then one can fit data to different parameters r characterizing 
the space-time structure of the source. 

One should realize at this point that there were already attempts to study the 
BEC with power-like (Lorenzian type) shape of source function: p(£) = 47I i R4 - — * 5//2 
which leads to C 2 = l + exp(-2RQ) (with £ 2 = (x 2 + y 2 + z 2 ) + (ct) 2 and Q defined 
as in our case) P0| . Such a form is nearest to our p{r) and, as it turns out, gives 
the best fit (in terms of the x 2_va hies) to data considered in [20|. However, differ- 



ences between this fit and other more conventional ones (i.e., based on gaussian or 
exponential shapes of the source) were not dramatic. This means, that in reality it 
will be very difficult to establish by means of BEC the possible existence of fractal 
structure of the emitting source. Perhaps the event-by-event analysis of data with 
some preselection of the initial conditions (in terms of energy, centrality, multiplicity 
etc.) will be necessary in order to perform such investigations. 

Our approach must be regarded as preliminary because of our choice of treatment 
of BEC. Notwithstanding its obvious deficiencies (already mentioned in [[!]]) it seems, 
however, fully adequate for the present study which is, as mentioned above, of only 
limited scope. However, even in such form it seems to indicate the relevance of the 
fact of a possible fractal structure of the space-time of the emitting region preventing 
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particles from different branches to be in the same emitting cell irrespectively of the 
smallness of differences in their positions or momenta. This stresses the problem of 



the distance in the cascade and the like, recently discussed in pi 
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Figure Captions: 



Fig. 1 The scheme of our cascade process. 

Fig. 2 Example of rapidity distributions of secondaries for totally symmetric (a) and 
totally asymmetric one-dimensional cascades with asymmetric (b) and sym- 
metric (c) emission of particles calculated for M = 40 GeV, jjl = 0.3. His- 
tograms are for fixed k: P(k) = S(k — 0.25). Open symbols display results 
for cascades with k\^ distributed randomly according to P(k) = (1 — k) a with 
a = 1 (in both cases multiplicities are the same: iV( ) ~ 11.5 and N^,c) — 4.5). 
Full lines present the most probable (one dimensional) thermal-like distribu- 
tions given by §) with (3 = 0.028 for (a) and /3 = -0.133 for (6) and (c), 
respectively, (calculated as in |T3| for iV( Q ) = 11.5 and iV(&, c ) — 4.5). 



Fig. 3 Multiplicity distributions P{N) for (one dimensional) cascades of masses M = 
10, 40, 100 GeV (for \x = 0.3 GeV and kip given by the same triangle dis- 
tributions P(k) as in Fig. 2): (a) - symmetric case with respective mean 
values (N) = 7.39, 14.67, 22.76 and dispersions a = 2.37, 5.12, 8.47; (&) - 
asymmetric case with (N) = 3.97,4.88, 5.48 and a = 1.07, 1.27, 1.41. 

Fig. 4 Example of 2 nd and 3 th scaled "horizontal" moments F\ as function of the 
number of bins nu n — Y jby (where Y is the rapidity range considered and Sy 
the bin size) for one-dimensional cascade of M = 40 GeV and \x = 0.3 GeV 
with k\ = k 2 = 0.25 and k 12 chosen randomly in the same way as in Fig. 3. 
The results for M = 100 GeV are essentially identical. 

Fig. 5 Density distribution of the production points p(r) for one- dimensional cascades 
(r = V* , fJ> = 0.3 GeV) - left panels, and for three-dimensional cascades 
(r = \/x 2 + y 2 + z 2 , /i = 0.14 GeV) - right panels. Two different choices of 
the evolution parameter r are considered: r = 0.2 fm - upper panels, and 
r = 0.2/M (in fm, the mass M is the parent mass in a given vertex) - lower 
panels. Each panel shows results for three different masses M of the source: 
M = 10, 40 and 100 GeV. In all cases k is chosen from the same triangle 
distribution distribution as in Fig. 2. 

Fig. 6 The C 2 (Q = \pi — Pj\) for the sources presented in Fig. 5: left panels - one- 
dimensional cascades, right panels - three-dimensional cascades; upper panels 

- time evolution parameter is set constant and equal r = 0.2 fm, lower panels 

- time evolution parameter is chosen as r = 0.2/M (in fm, the mass M is the 

13 



parent mass in a given vertex). Each panel shows results for three different 
masses M of the source: M = 10, 40 and 100 GeV. In all cases k is chosen 
from the same triangle distribution distribution as in Fig. 2. 
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Fig.5 
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Fig.6 



